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Using the anti-de Sitter/conformal field theory correspondence, we relate the shear viscosity n of 
the finite-temperature M = 4 supersymmetric Yang-Mills theory in the large N, strong-coupling 
regime with the absorption cross section of low-energy gravitons by a near-extremal black three- 
brane. We show that in the limit of zero frequency this cross section coincides with the area of the 
horizon. From this result we find r\ = jN 2 T 3 . We conjecture that for finite 't Hooft coupling <?Y M iV 
the shear viscosity is r\ = f(g\ M N)N 2 T 3 , where f(x) is a monotonic function that decreases from 
0(x~ 2 ln _1 (l/x)) at small x to 7r/8 when x — > oo. 



Introduction. — At finite temperatures, the large dis- 
tance, long time behavior of gauge theories is described, 
as in any other fluid, by a hydrodynamic theory Q. To 
write down the hydrodynamic equations one has to know 
the thermodynamics (i.e., the equation of state) of the 
medium, as well as the transport coefficients: the shear 
and the bulk viscosities, the electrical conductivity [in 
the presence of a U(l) gauge group], and the diffusion 
constants (in the presence of conserved global charges). 
Knowledge of these quantities in hot gauge theories is 
crucial for numerous applications, the most notable of 
which belong to the physics of the electroweak phase 
transition in the early Universe and of the quark-gluon 
plasma possibly created in heavy-ion collisions || . 

When the gauge coupling is small (which requires, in 
the case of QCD, temperatures much larger than the con- 
finement scale), both the equation of state and the trans- 
port coefficients are calculable perturbatively. At strong 
coupling (i.e., at temperatures not much larger than the 
confinement scale), thermodynamics can be found non- 
perturbatively by lattice simulations, but transport co- 
efficients are beyond the reach of all current numerical 
techniques. This situation is very unfortunate, since the 
quark-gluon plasma one hopes to create in heavy-ion ex- 
periments has relatively low temperature at which the 
perturbation theory works very poorly. 

Lacking methods to reliably compute the transport co- 
efficients of finite-temperature QCD, one should try to 
gain insight into models where these coefficients can be 
computed nonperturbatively. Recently, powerful tech- 
niques based on the anti-de Sitter/conformal field theory 
(AdS/CFT) correspondence have been developed, estab- 
lishing, in particular, the connection between the JV = 4 
supersymmetric Yang-Mills (SYM) theory in the large 
coupling, large N limit and classical ten-dimensional 
gravity on the background of black three-branes j| 
This allows one to perform analytical calculations in a 
strongly coupled four-dimensional gauge theory. 

In this Letter, we compute the shear viscosity 77 of the 



strongly coupled finite-temperature JV = 4 SYM theory 
(the bulk viscosity of this theory vanishes due to scale in- 
variance). We first relate, using previously known results 
from the AdS/CFT correspondence, the shear viscosity 
with the absorption cross section of low-energy gravitons 
falling perpendicularly onto near-extremal black three- 
branes. We further show that this cross section is equal 
to the the area of the horizon, in a way very similar to 
the case of black holes || . These facts provide enough in- 
formation for us to find that 77 = ^N 2 T 3 , provided both 
the 't Hooft coupling and N are large. Remarkably, the 
shear viscosity approaches a constant value in the large 
't Hooft coupling limit, and this value is related to the 
area of the horizon of the black brane. 

The viscosity. — To start our discussion, we briefly re- 
view the notion of viscosity in the context of finite- 
temperature field theory. Consider a plasma slightly out 
of equilibrium, so that there is local thermal equilibrium 
everywhere but the temperature and the mean velocity 
slowly vary in space. We define, at any point, the local 
rest frame as the one where the three-momentum den- 
sity vanishes: T oi = 0. The stress tensor, in this frame, 
is given by the constitutive relation, 

Tij = 5 tj p -rjfdiUj + d 3 Ui - ^StjdkU^J 

-(SijdkUk , (I) 

where Ui is the flow velocity, p is the pressure, and 77 
and £ are, by definition, the shear and bulk viscosities 
respectively. In conformal field theories like the JV = 4 
SYM theory, the energy momentum tensor is traceless, 
T M M = 0, so s = Too = 3p an d the bulk viscosity vanishes 
identically, £ = 0. 

All kinetic coefficients can be expressed, through Kubo 
relations, as the correlation functions of the correspond- 
ing currents || . For the shear viscosity, the correlator is 
that of the stress tensor, 

77= lim-^ [dtdxe iut ([T xy (t,x.),T xy (0,0)]) 
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lim— [GaH-GrH], 



(2) 



where the average (. . .) is taken in the equilibrium ther- 
mal ensemble, and Ga and Gr, are the advanced and 
retarded Green functions of T xy , respectively. In Eq. (|^), 
the Green functions are computed at zero spatial mo- 
mentum. Though Eq. (g) can, in principle, be used to 
compute the viscosity in weakly coupled field theories, 
this direct method is usually very cumbersome, since it 
requires rcsummation of an infinite series of Feynman 
graphs. This calculation has been explicitly carried out 
only for scalar theories (lcf] , A more practical method is 
to use the kinetic Boltzmann equation, which gives the 
same results as the diagrammatic approach |fll| . 

For gauge theories at weak coupling, g 2 N <C 1, where 
throughout this paper g = <?ym is the gauge coupling, the 
shear viscosity has the following parametric behavior, 



v = c 



N 2 T 3 



(g 2 N) 2 \n(l/g 2 N) 



(3) 



Basically, ij is proportional to the product of the energy 
density e ~ N 2 T A and the transport mean free time r ~ 
[{g 2 N) 2 \n(l/g 2 N)T}- 1 . The numerical coefficient G in 
Eq. (^), in principle, can be computed by solving the 
linearized Boltzmann equation |l2| . 

The relation to graviton absorption. — The key obser- 
vation underlying this work is that the right hand side 
of the Kubo formula (||) is known to be proportional to 
the classical absorption cross section of gravitons by black 
three-branes |I^,[l4|]. For completeness, we recall here the 
basic argument leading to this correspondence. Consider, 
in type IIB string theory, a configuration of N D3-branes 
stacked on top of each other. The low-energy theory liv- 
ing on the branes is the Af = 4 U(JV) SYM theory. On 
the other hand, if N is large, the stack of D3-branes 
has large tension, which curves space-time. In the limit 
of large 't Hooft coupling g 2 N , the three-brane geometry 
has small curvature and can be described by supergravity. 
Therefore, we have two descriptions of the same physics 
in terms of strongly coupled gauge theory on the branes 
and classical gravity on a certain background. 

If one sends a graviton to the brane, there is some prob- 
ability that it will be absorbed. On the gravity side, the 
absorption cross section can be calculated by solving the 
wave equation on the background metric. On the gauge 
theory side, the rate of graviton absorption measures the 
imaginary part of the stress tensor-stress tensor corre- 
lator, since gravitons polarized parallel to the brane are 
coupled to the stress-energy tensor of the degrees of free- 
dom on the brane. The relation between the absorption 
cross section o~(u>) of a graviton with energy u>, polar- 
ized parallel to the brane (say, along the xy directions) 
and falling at a right angle on the brane is related to the 
correlator in field theory as Pq,pT[ 



tr(w) - — / dtdxe iut ([T xy (t,x),T xy {0,0)}) , (4) 



where k = Vo7rG, G being the ten-dimensional gravita- 
tional constant. The relation ([!]) has been explicitly ver- 
ified for zero-temperature field theory (or extremal black 
branes, in the gravity language) |13ll4l. Such a check is 
possible because there is a nonrenormalization theorem 
for the correlator of the stress-energy tensor Jl9]. 

At finite temperature T, Eq. (|j) relates the graviton 
absorption cross section by a near-extremal black brane 
having the Hawking temperature equal to T with a corre- 
lator in the hot SYM theory [|l6| . Since there is no super- 
symmetry and no nonrenormalization theorem is known 
to work at finite temperature, one cannot explicitly verify 
the relation (|J). We instead view Eq. (f|) as a prediction 
of theory. In particular, taking the to — > limit one 
can relate a(u> — 0) with the shear viscosity of hot SYM 
plasma, 



1 = ^"(0) 



(5) 



Equation ([5|) implies that, for nonextremal black branes, 
the graviton absorption cross section must not vanish in 
the limit of zero frequency (in contrast to the extremal 
case where o-(ui) ~ ui 3 at small ui jl3 14|), and, by com- 
puting the zero- frequency value of a one obtains the shear 
viscosity of the hot SYM plasma. The problem of com- 
puting the shear viscosity is now reduced to a problem 
of classical gravity. 

The metric of a nonextremal black three-brane has the 

form Ir^Jrl 



ds 2 = H- 1/2 (r)[-f{r)dt 2 + dx 
+H 1/2 {r)[f- 1 {r)dr 2 + r 



2 ' ™ 2 dn 2 ] 



(0) 



where H(r) = 1 + i? 4 /r 4 and f(r) = 1 - r 4 /r 4 . The 
extremal case corresponds to ro = 0; the limit relevant 
for us is the near-extremal one, ro <C R. This metric has 
a horizon at r = ro. From the existence of this horizon 
one should expect cr(0) to be nonvanishing. Running 
ahead, we will show, by solving the wave equation on the 
metric (||), that er(0) is equal to the area of the horizon, 



cr(0) = TrVgi? 2 



(7) 



(the numerical coefficient ir 3 is simply the area of the 
unit five-sphere). Using the formula for the Hawking 
temperature of the metric (Jq), 



T 



ttR 2 ' 



(8) 



and the relation between R and N which is obtained 
by identifying the Arnowitt-Deser-Misner mass per unit 
volume of the three-brane with the tension of a stack of 
N D3-branes [ftl, 
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i? 4 = 



kN 



2tt 5 /2 ' 

we find the shear viscosity to be 



V 



2 rp3 



-NT 



(9) 



(10) 



This is the main result of the paper. Up to a constant, 
the shear viscosity is equal to the entropy density Jl9],^0| . 
Both quantities are proportional to the area of the hori- 
zon. 

Solution to the radial wave equation. — Now let us 
show how Eq. ^ is obtained. We have to solve the 
s-wave radial equation for a minimally coupled scalar 
(such as the graviton polarized parallel to the brane), 
dtii^f—gg^d^cf) = 0. In the metric (||) this equation 
acquires the form 



5r 4 



r(r 4 



( r 4 



A\2 



4> = o. 



(ii) 



The method we use to solve Eq. ( |ll[ ) is the same matching 
method that was used in the extremal case jl3|,[l4| . Since 
ultimately we are interested in the limit lo — > 0, we will 
assume ui <C T . More details about our method, as well 
as the solution to the radial equation in the opposite limit 
w»T and for higher partial waves, can be found in Ref. 
pl[ . Earlier attempts to compute the absorption rate 
by nonextremal black branes were made in Refs. p2[ . 
As in the extremal case, we search for the solution in 
several regions and match the result wherever the regions 
overlap. Let us go from small r to large r, starting from 
the horizon r = tq. The first region is the one just outside 
the horizon: r > ro, r — ro <C vq. In this case Eq. (|I7 
has the form 



A 2 



'"o 



16 (r — ro)' 



= 0. 



(12) 



where A = w/(irT) < 1. The solution to Eq. (|lg) is 

-iA/4 



= A 1- 



(13) 



where we have chosen the sign of the exponent so that the 
solution corresponds to an incoming wave at the horizon. 
When A is small, (|lj) is basically a constant, A, except 
for an exponentially small region near ro. In the next 
region, ro < r <C lo^ 1 (excluding r exponentially closed 
to ro), the term proportional to uj 2 in the left-hand side 
of Eq. (Thy) can be dropped. Indeed, 



cA 8 (r 4 -roT^«(r-r )- 



(14) 



due to r -C u> , and 



w Vi? 4 (r 4 - r^)- 2 « (r - ro)" 2 (15) 
since lj -C T ~ r i?~ 2 . Equation (0) now has the form 



5r 4 



'o 



r(r 4 — Tq) 
which possesses a trivial solution, 
J) = A, 



' = 0. 



(16) 



(17) 



which matches smoothly with Eq. (|l3|). Finally, in the 
outermost region, r 3> R 3> ro, Eq. (pd|) is simplified to 



d 2 4> 5dcj) 2 , 
dr z r dr 



(18) 



which can be solved in terms of the Bessel functions, 

^)=a#^+/3^M, r»R. (19) 



(wr) 



(wr) 



The regimes of validity of Eq. (|17j), ro < r <C and of 
Eq. (|19|), r>ij, has an overlap since w _1 ^> i? (this is 
the consequence of R 3> ro and w <C T) . In order for Eq. 
© to match with Eq. (0) 

in the overlapping region, 

one should require 



= 8A 



(3 = 0. 



(20) 



The field at large distances can be decomposed into an 
incoming wave and an outgoing wave, 



-) = AA 



(wr) 2 (cor) 2 



(2)/ 



(21) 



The absorption probability P is the ratio of the flux 
at r = r from Eq. ( |l3| ) and the flux from the incoming 
wave m Eq. @). We find 



7T 

32" 



P = —uj^rlR 2 . 



(22) 



Since the absorption cross section a is related to P by 



32^ 2 



(23) 



we arrive to Eq. (0) , which coincides with the area of the 
horizon. This is very similar to the universal result for 
black holes @. 

Notice that in deriving Eq. (ffl) we require to be 
much smaller than the Hawking temperature. The ab- 
sorption cross section will deviate substantially from the 
zero-frequency limit if ui is of order T. In particular, 
the next correction to Eq. (]?]) is of order uj 2 /T 2 with a 
computable coefficient pl[ . 

Discussion. — We have shown that the shear viscosity 
can be computed in the strongly coupled Af = 4 SYM 
theory from the AdS/CFT correspondence. Now let us 
try to interpret the result (|Toj) . The power of T in r/ is 
completely fixed by the dimensionality of r/ and the scale 
invariance of the theory. The factor N 2 apparently comes 
from the number of degrees of freedom in the plasma. It 



3 



is remarkable that the shear viscosity approaches a con- 
stant value as one sends the 't Hooft coupling to infinity. 
From the relation r) ~ er, one can interpret this behavior 
as the indication that the "relaxation time" r remains of 
order T _1 (but not much smaller) in the strong coupling 
limit. Since the inverse relaxation time is comparable 
to the energy per degree of freedom, the strongly cou- 
pled plasma cannot be viewed as a collection of particles, 
and the formula 77 ~ er does not applies in the strict 
sense. However, one should expect that the counting of 
the powers of N still works. This counting, combined 
with the expressions for 77 is the weak-coupling [Eq. (^)] 
and strong-coupling [Eq. (0)] limits, suggests that for 
finite 't Hooft coupling g 2 N the shear viscosity has the 
form 



77 = f(g 2 N)N 2 T 3 



(24) 



where f(x) ~ x~ 2 In - (1/a;) when x <C 1 and f(x) = tt/8 
when x ^> 1. It is most likely that f(x) is a monotonic 
function of x. One way to verify this conjecture is to com- 
pute the 0(l/g 2 N) correction to 77 in the strong coupling 
limit. If f(x) is monotonic, then this correction must be 
positive. This is analogous to the behavior of the free 
energy po| , except that at small coupling 77 — > 00, while 
the free energy remains finite. 

Recalling that cr(ui) deviates substantially from <r(0) 
when u> ~ T, we see that the hydrodynamic theory 
can describe processes occuring during times much larger 
than T _1 , but breaks down for those whose typical time 
scale is of order or less than T -1 . One also should ex- 
pect hydrodynamics to work at spatial distances much 
larger than T" 1 , but not at distances of order or less 
than T _1 . This is consistent with T _1 playing the role 



of the relaxation time in the limit g N 



There is 



apparently no separation of scales in the strong coupling 
regime that would make a kinetic description possible: 
T _1 is the only time/length scale. Thus, the viscosity 77 
cannot be computed from a Boltzmann-type equation. 

In this paper, we have confined our attention to the 
most important transport coefficient — the shear vis- 
cosity. As mentioned above, the bulk viscosity vanishes 
identically due to the exact scale invariance of the Af = 4 
SYM theory. It would be useful to compute other trans- 
port coefficients in this theory (for example, the diffusion 
constant of the R charges) at strong coupling using the 
AdS/CFT correspondence. 
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